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The  compression  and  fatigue  behavior  of  a  Ca65Mg15Zn2o  bulk-metallic  glass  (BMG)  was 
studied  in  air  at  room  temperature.  During  the  preparation  of  cubical  samples  of  the 
Ca65Mg15Zn2o  for  compression  and  fatigue  investigations,  small  spherical  cavities  were  found. 
Under  both  monotonic  and  cyclic  compression  loadings  of  the  samples,  fractures  initiated  at 
these  cavities  and  propagated  in  a  direction  generally  parallel  to  the  loading  axis.  Finite-element 
analysis  (FEA)  was  used  to  model  the  fracture  behavior.  The  FEA  of  a  centrally  located 
spherical  void  showed  that  under  compression  loading,  large  tensile  stresses  evolved  in  the 
cavities.  The  orientation  of  the  maximum  principal  stress  (PI)  was  found  to  be  normal  to  the 
direction  of  crack  propagation,  which  is  consistent  with  the  experimental  finding.  Stresses  in 
deeply  embedded  adjacent  voids  and  those  in  superficial  voids  were  also  studied.  The  influence 
of  the  void  location  in  the  cubical  sample  on  the  fracture  behavior  was  quantitatively  discussed. 
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I.  INTRODUCTION  AND  OUTLINE 
OF  THE  PROBLEM 

Bulk-metallic  glasses  (bmgs)  lack  grain 
boundaries  and  any  long-range  order.  Therefore,  they 
are  significantly  different  from  crystalline  materials.  The 
BMGs  have  many  excellent  physical  and  chemical 
properties,  such  as  high  strengths,  high  elastic  limits, 
high  wear  and  corrosion  resistances,  almost  perfect 
as-cast  surfaces,  and  good  fracture  toughness. [1,21  Since 
BMGs  can  be  applied  as  potential  structural  materials, 
the  mechanical  behavior  of  BMGs  is  being  studied 
extensively. [3  8] 

The  Ca-based  BMGs  have  many  unique  properties, 
such  as  low  density  (~2.0  g/cc),  low  Young’s  modulus 
(~20  to  30  GPa),  low  shear  modulus  (~8  to  15  GPa),  low 
glass-transition  temperature  (Tg  ~  100  °C  to  190  °C), 
and  low  crystallization  temperature  ( Tx  ~  130  °C  to 
240  °C).[8  Moreover,  Ca-based  BMGs  have  very 
good  glass-forming  ability  (GFA).  The  elastic  modulus 
of  Ca-based  BMGs  is  comparable  to  that  of  human 
bones,  and  Ca,  Mg,  and  Zn  are  biocompatible. [11]  These 
features  make  the  Ca-Mg-Zn-based  alloys  attractive  for 
use  in  biomedical  applications. [10]  Since  Ca-based 
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BMGs  were  successfully  synthesized^12,1 3]  Ca-based 
BMGs  are  attracting  more  and  more  attention. 

In  a  previous  article, [14]  the  monotonic  compression 
and  compression-compression  fatigue  experiments  were 
employed  to  investigate  the  mechanical  properties  of  the 
Ca65Mgi5Zn2o  (at.  pet)  BMG  at  room  temperature  in 
air.  Cubical  samples,  4  mm  on  a  side,  were  prepared. 
Under  compression  loading,  it  was  observed  that 
fractures  originated  from  small  deeply  embedded  spher¬ 
ical  voids,  typically  25  to  50  /an  in  diameter,  which 
functioned  as  crack-nucleation  sites.  The  alloy  was 
macroscopically  brittle  and  crack  propagation  occurred 
in  a  direction  nearly  parallel  to  the  loading  axis,  which 
demonstrated  a  dominant  splitting  failure  with  fracture 
planes  from  0  to  20  deg  relative  to  the  loading  axis. 
Figure  1  demonstrated  an  example  of  splitting  fracture 
under  monotonic-compression  loading.  In  general,  this 
splitting  fracture  is  observed  in  brittle  materials,  includ¬ 
ing  Fe-based  BMGs.[14]  This  phenomenon  could  be 
caused  by  the  brittleness  of  materials  and  voids  in  the 
materials.  To  explain  the  observed  splitting  behavior, 
finite-element  analysis  (FEA)  was  conducted  in  the 
current  work. 


II.  FINITE-ELEMENT  MODELING 

In  the  finite-element  method, [15]  a  structure  is  subdi¬ 
vided  into  a  number  of  finite-sized  volumes,  called 
elements.  In  each  element,  the  displacements,  strains, 
and  stresses  are  approximated  by  polynomial  shape 
functions.  Connections  between  elements  are  made 
through  a  finite  number  of  grid  points,  commonly 
called  nodes.  Interactions  among  elements  occur  solely 
through  the  forces  exerted  at  the  grid  points.  Element 
material  properties  and  a  structure’s  geometry  are  used 
to  develop  the  structure’s  stiffness  characteristics. 
Known  loads  are  applied  at  nodes  to  the  structure 
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Fig.  1  SEM  fractographs  demonstrating  (a)  a  porelike  flaw,  microcrack, 
resulted  in  macroscopic  splitting  failure  under  compression. [14] 


with  a  given  stiffness  and  the  strain  and  stress  distri¬ 
butions  in  the  structure  are  calculated.  The  develop¬ 
ment  of  the  finite-element  model  was  straightforward. 
The  necessary  elements  are  the  geometry,  material 
properties,  boundary  conditions,  loading,  and  element 
formulation. 

The  first  case  studied  was  that  of  a  deeply  embedded 
(2  mm  from  each  face)  centrally  located  void  of  50- fim 
diameter,  which  shows  both  the  solid  model  and  the 
finite-element  model.  Results  of  this  analysis  demon¬ 
strate  the  essential  nature  of  the  stress  distribution  in  the 
void.  In  addition,  different  failure  theories  could  be 
applied  to  determine  which  would  best  explain  the 
observed  behavior  of  dominant  splitting  fracture.  It  was 
anticipated  that  the  results  of  this  analysis  would,  in 
addition  to  determining  general  stress  levels,  provide 
some  insight  into  the  void  behavior,  relative  to  the  size 
and  shape  of  the  cubical  sample. 

Case  2  was  the  analysis  of  three  inline  50-^m-diameter 
voids,  with  a  center-to-center  spacing  of  150  fim.  The 
results  of  this  analysis  allowed  the  effect  of  one  void 
upon  another  to  be  assessed.  It  would  show  whether  the 
stress  concentration  at  one  void  had  a  meaningful  effect 
upon  the  stress  state  in  a  nearby  adjacent  void. 

Case  3  examined  the  stress  distribution  in  a  hemi- 
spherically  shaped  surface  void.  Although  no  fractures 
emanating  from  these  voids  were  observed,  it  was 
nevertheless  deemed  important  to  learn  if  the  stress 
distribution  in  superficial  voids  was  fundamentally 
different  from  that  of  a  deeply  embedded  void  or  if 
surface  voids  posed  different  risks. 

The  fourth  and  last  case  permitted  an  examination  of 
the  size  effect.  It  was  identical  to  case  1  with  the 
exception  that  the  void  diameter  was  100  fi m,  rather 
than  50  fim.  If  the  stress  state  in  a  100-/im-diameter  void 
was  to  be  more  severe  than  that  in  a  50-^m-diameter 
void,  then  compression-compression  fatigue  perfor¬ 
mance  could  be  improved  by  controlling  the  void  size. 

The  FEA  model  used  in  case  1  was  a  cube,  2  mm  on  a 
side,  with  one  octant  of  the  centrally  located  void 
remaining.  Since  there  are  three  planes  of  the  geometrical 
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and  loading  symmetry,  it  was  only  necessary  to  model 
one  octant. 

The  Ca65Mg15Zn2o  is  quite  brittle  and  exhibits  a 
purely  linear  stress/strain  response  up  to  the  loading 
strains  of  ~1.5  pet.  It  has  a  Young’s  modulus  of 
20  GPa,[14]  fracture  strength  of  364  MPa,[14]  and 
Poisson’s  ratio  of  0.30.[16] 

Applied  boundary  conditions  were  zero  normal  dis¬ 
placements  at  all  nodes  on  each  plane  of  the  symmetry. 
Model  loading  in  all  cases  was  a  compressive  displace¬ 
ment  of  0.0254  mm  applied  in  the  x  direction.  This 
arrangement  is  equivalent  to  compressive  strain  of 
1.27  pet. 

In  the  case  2  finite-element  model,  due  to  the 
symmetry  in  both  the  geometry  and  loading,  only  one 
octant  is  required  for  the  solution,  since  no  displace¬ 
ments  normal  to  the  surfaces  that  define  the  symmetry 
can  occur.  Thus,  only  one-and-one-half  voids  are 
modeled. 

Isoparametric  brick  elements  with  quadratic-shape 
functions  were  used  throughout.  NEi  Nastran,[15]  a 
displacement-based  general  purpose  finite-element  code, 
well  suited  to  problems  of  this  type,  was  used  for  the 
solution. 


III.  RESULTS 

A  partially-zoomed  view  of  the  case  1  maximum 
principal  stress  distribution  is  shown  in  Figure  2.  The 
stress  peak  value  is  180.28  MPa.  Its  location  is  indicated 
as  point  A  in  Figure  10.  Point  B,  which  is  located  on  the 
surface  of  the  void  along  the  z-axis  and  has  a  zero  stress 
level,  is  also  noted  in  Figure  2.  The  maximum  principal 
stress  (PI)  distribution  along  arc  A-B  is  shown  in 
Figure  3.  In  Figure  4,  the  same  stress  distribution  is 
presented,  but  the  viewing  direction  is  parallel  to  the 
loading  plane.  The  singular  feature  seen  in  Figure  4  is 
the  circularity  of  the  stress  distribution.  It  is  axisym- 
metric  with  respect  to  the  x-axis.  The  deformed  shape 
of  the  void  is  generally  ellipsoidal,  in  which  the 
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Fig.  2  Zoomed  view  of  the  maximum  principal  stress  distribution 
in  a  deeply  embedded  50  /an  diameter  spherical  void  (case  1).  Arc 
A  B  has  been  identified.  Units  are  MPa. 


Fig.  4  Maximum  principal  stress  distribution  in  a  deeply  embedded 
50  /nn  diameter  spherical  void  (case  1).  The  viewing  direction  is 
parallel  to  the  applied  load.  Units  are  MPa. 
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Fig.  3  Graph  of  the  maximum  principal  stress  distribution  along 
arc  A  B  in  a  deeply  embedded  50  /an  diameter  spherical  void 
(case  1).  Figure  2  provides  a  definition  of  arc  A  B. 


Fig.  5  Zoomed  view  of  the  maximum  principal  stress  distribution 
in  case  2.  Note  a  similarity  to  case  1  results.  Also  note  that  the  void 
shapes  are  different  because  there  are  three  planes  of  symmetry  in 
the  centrally  located  void,  but  only  two  planes  of  symmetry  in  the 
void  whose  center  is  at  an  x  coordinate  of  150  /an.  Units  are  MPa. 


deformation  is  superimposed  over  the  original  geome¬ 
try,  which  is  shown  in  bold  lines. 

Case  2  results  are  shown  in  Figures  5  and  6.  Figure  5 
is  a  zoomed  view  of  the  maximum  principal  stress 
distribution  in  the  vicinity  of  two  adjacent  voids. 
Figure  6  is  a  plot  of  the  maximum  principal-stress 
distribution  along  arc  A-B-C.  The  maximum  PI  values 
were  observed  to  be  176.11  MPa,  which  is  the  same  as 
case  1.  Only  small  void-to-void  differences  (6.93  MPa) 
were  found. 

Case  3  results  are  seen  in  Figure  7.  The  maximum  PI 
value  was  160.69  MPa,  which  is  smaller  than  the  value 
observed  in  case  1.  In  surface  voids,  the  maximum 
principal  stress  distribution  is  not  axisymmetric,  but 
shows  two  pronounced  stress  concentrations.  A  graph 
of  the  maximum  principal  stress  distribution  passing 
through  the  two  stress  concentrations  is  presented  in 
Figure  8. 
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Maximum  Principal  Stress  Along  Arc  A-B-C 


Angular  Position  (A=0  deg,  C=180  deg) 

Fig.  6  Maximum  principal  stress  distribution  along  arc  ABC 
(case  2). 
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50-Micron-Diameter  vs  100-Micron-Diameter 
Maximum  Principal  Stresses 


Angular  Position  (A  -  0  deg,  B  =  90  deg) 


Fig.  7  Zoomed  view  of  the  maximum  principal  stress  distribution 
for  case  3.  Note  that  the  stress  distribution  does  not  manifest  the 
kind  of  symmetry  seen  in  other  cases  and  that  the  peak  value  is  less 
than  that  seen  in  other  cases.  Units  are  MPa. 


Fig.  9  Comparison  of  the  maximum  principal  stress  distributions  in 
cases  1  and  4  along  arc  A  B,  which  is  defined  in  Fig.  2.  Note  that 
the  values  are  very  nearly  identical  at  every  nodal  point.  This  fact 
confirms  the  conclusion  that  for  small  voids,  PI  values  are  indepen 
dent  of  the  void  diameter. 


Maximum  Principal  Stress  Distribution  In 
50  Micron  Diameter  Hemispherical  Surface 
Void  Along  Arc  A-E3-C 


Fig.  8  Graph  of  the  maximum  principal  stress  distribution  along 
arc  ABC  (case  3).  Note  the  symmetry  in  the  distribution  and  the 
fact  that  peak  value  is  less  than  that  in  other  cases. 

Case  4  results  are  identical  to  those  of  case  1.  Figure  9 
is  a  graph  that  compares  the  PI  stress  distribution  in  a 
50-^m-diameter  void  to  that  in  a  100-/im-diameter  void. 
Since  the  FEA  mesh  in  both  cases  contains  the  same 
number  of  elements  in  the  void,  the  results  reflect  PI 
values  at  corresponding  angular  locations  along  arc 
A-B,  which  is  defined  in  Figure  2.  There  are  essentially 
identical  results  for  the  50-  and  100-/im  voids  in 
Figure  9. 


IV.  DISCUSSION 

From  the  case  1  results  (Figures  2  through  4),  we 
learn  that  the  PI  distribution  is  axisymmetric.  Since  the 
specimen  is  a  cube,  we  can  conclude  that  the  PI 
distribution  is  unaffected  by  the  general  part  shape. 


This  trend  is  not  surprising,  given  that  the  overall 
dimensions  of  the  part  are  so  much  larger  than  the  void. 
Thus,  the  void  behaves  as  though  it  were  situated  in  a 
semi-infinite  medium.  The  deformed  geometry  is  gener¬ 
ally  ellipsoidal,  with  equal  positive  displacements  occur¬ 
ring  along  the  y-  and  z-axes,  and  a  negative 
displacement  along  the  z-axis.  This  change  in  the  shape 
leads  to  a  flattening  of  the  void  in  the  vicinity  of  the 
loading  axis,  resulting  in  localized  stretching  that  causes 
large  transverse  tensile  stresses.  It  is  this  tensile  stress 
that  is  responsible  for  the  observed  splitting  behavior. 

Since  the  material’s  purely  linear  stress-strain 
response  suggests  a  brittle  failure  mechanism,  it  is 
natural  to  apply  the  Rankine,  or  the  maximum  normal 
stress  theory  of  failure. [1 7,1 8]  Indeed,  fractures  do  initiate 
at  or  very  close  to  the  point  of  the  maximum  principal 
stress,  as  predicted  by  the  model.  Any  differences 
between  the  observed  behavior  and  that  predicted  by 
the  model  are  likely  attributable  to  small  deviations 
from  the  sphericity  in  actual  voids.  Crack  propagation 
occurs  in  a  direction  parallel  to  the  loading,  i.e.,  along 
the  positive  x-axis.  The  principal  direction  associated 
with  PI  is  within  5.83  deg  of  being  normal  to  the  x-axis. 
Thus,  the  calculated  initial  fracture  plane  should  be 
nearly  parallel  to  the  x-axis,  which  is  consistent  with  the 
experimental  observation  (Figure  1).[14] 

From  the  case  2  results  (Figures  5  and  6),  it  is  seen 
that  voids  less  than  50  fix n  in  diameter  that  are  separated 
from  each  other  by  more  than  150  /an  are  unlikely  to 
cause  increases  in  the  stress  above  levels  found  in 
isolated  voids.  PI  stress  levels  are  equivalent  to  those  of 
case  1.  The  principal  direction  associated  with  PI  is 
within  1.46  deg  of  being  normal  to  the  x-axis.  The 
calculated  initial  fracture  plane  should  be  almost  par¬ 
allel  to  the  x-axis,  as  observed  experimentally  in 
Figure  1. 

The  case  3  results  (Figures  7  and  8)  suggest  that 
fractures  should  preferentially  initiate  in  deeply  embed¬ 
ded  voids,  rather  than  in  superficial  voids,  since  the 
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maximum  PI  values  are  lower  in  surface  voids  than  in 
deeply  embedded  voids.  In  fact,  no  fractures  have  been 
observed  to  have  nucleated  at  surface  voids.  The 
principal  direction  associated  with  PI  is  within 
0.21  deg  of  being  normal  to  the  x-axis.  Thus,  the 
calculated  initial  fracture  plane  should  be  almost  par¬ 
allel  to  the  x-axis,  which  is  consistent  with  the  experi¬ 
mental  results  (Figure  1,  Reference  14). 

The  case  4  stress  distribution  (Figure  9),  although 
identical  to  that  of  case  1,  is  nevertheless  important.  For 
voids  less  than  100  /im  in  diameter  in  a  4-mm  cube,  the 
results  strongly  suggest  that  the  stress  distribution  is 
independent  of  the  diameter.  This  surprising  result  is 
consistent  with  the  case  1  conclusion  that  the  void  is 
located  in  a  semi-infinite  medium.  It  is  very  similar  to 
the  well-known  elasticity[19,20]  result  of  a  circular  hole  in 
an  infinite  plate  under  plane  stress  uniaxial  loading.  In 
the  elasticity,  stress  results  on  the  surface  of  the  circular 
void  are  independent  of  the  size  of  the  void  and 
maximum  principal  stress  in  response  to  compression 
loading  occurs  along  the  load  line.  The  implication  of 
this  conclusion  on  any  effort  to  improve  fatigue  perfor¬ 
mance  is  that  all  voids,  regardless  of  size,  should  be 
eliminated.  It  is,  thus,  insufficient  merely  to  reduce  the 
void  size  because  the  stress  state  cannot  be  changed. 


V.  SUMMARY 

This  study  focused  on  explaining  the  splitting  behav¬ 
ior  of  Ca  BMGs  under  compression  loading.  To 
accomplish  this  task,  finite-element  analyses  were  used 
to  determine  the  stress  states  in  small  spherical  voids, 
and  the  maximum  normal-stress  theory  of  failure  was 
applied.  The  FEA  results  were  used  to  accurately  predict 
the  crack-nucleation  sites  and  the  initial  direction  of 
crack  propagation,  which  is  generally  parallel  to  the 
loading  direction,  as  observed  experimentally.  It  was 
found  that  adjacent  voids  do  not  interact  with  each 
other,  and  cracks  should  not  nucleate  from  small  surface 
voids.  The  most  important  result  was  that  the  stress 
distribution  is  independent  of  the  void  diameter.  This 
trend  led  to  the  conclusion  that  improvements  in  the 
fatigue  performance  should  not  be  expected  by  reducing 
the  void  size. 
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